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Abstract. We present XMM- Newton observations of the 
bright quasar PKS 0558-504. The 0.2-10 keV spectrum 
is dominated by a large, variable soft X-ray excess. The 
fastest flux variations imply accretion onto a Kerr black 
hole. The XMM-Newton data suggest the presence of a 
'big blue bump' in PKS 0558—504 extending from the op- 
tical band to ~ 3 keV. The soft X-ray spectrum shows 
no evidence for significant absorption or emission-line fea- 
tures. The most likely explanation for the hot big blue 
bump is Comptonization by the multi-temperature corona 
of a thermal accretion disc running at a high accretion 
rate. 

Key words: galaxies: active - X-rays: galaxies - accre- 
tion discs - quasars: individual: PKS 0558—504 



1. Introduction 

The most prominent feature in the broad-band spectral 
energy distribution of many active galactic nuclei (AGN) 
is the 'big blue bump' (BBB), which contains a large frac- 
tion of the bolometric luminosity. In most broad-line AGN 
the BBB begins around 1 zxm, rises to a peak in the ultra- 
violet and then d ecline s again into the unobservable EUV 
(e.g. Elvis et al. 



1994). The soft X-ray excesses seen in 



many AGN may b e the high-energy tai l of th e BBB (e.g. 
Turn er & Pounds |1989| ; Walter et al. |1994| ; Laor et al. 
1997). The harder 2-10 keV X-ray spectrum can usually 



be quite well described by a power-law with photon-index 
T « 1.9. Most AGN also show an iron Ka fluorescence fea- 
ture around 6-7 keV and/or a hump of emission above 10 
keV (e.g. Pounds et al. 199C), widely interpreted in terms 
of Compton reflection of hard X-rays by optically-thick 
matter probably in the form of an accretion disc (George 
& Fabian |199l)). 



Three basic classes of model have been proposed to ex- 
plain the BBB. The first class involves thermal emission 
from an optically thick accretion disc surrounding the cen- 
tral black hole in which the primary heating is by viscous 
dissipation within the disc itself (e.g. Shields 1979| ). The 
disc thermal spectrum may be comptonized by a hot disc 
corona; a process which changes the observed shape of the 
continuum, helping to explain soft X-ray emission from 
discs which are otherwise too cool (e.g. Czerny & Elvis 
1987). A second class of models also involves thermal emis- 



Send offprint requests to: P.T. O'Brien 
Correspondence to: pto@star.le.ac.uk 



sion from an accretion disc, but with irradiation by a hard 
X-ray power-law (perhaps arising from magnetic flares) as 
the primary energy source (e.g. Ross & Fabian 1993). In 
these models the disc reprocesses the hard X-ray emission 
into thermal emission which is re-radiated at the local disc 
temperature. Irradiated discs may have flatter radial tem- 
perature profiles, and hence different spectra, compared 
to discs dominated by viscous dissipation. The final class 
of model involves optically thin, t herm al free-free emis- 
sion from a hot plasma (Barvainis 1992), although single 
temperature models of this kind do not fit the observed 
continuum very well (Siemiginowska et al. 1995). 

While the spectral form of the optical and UV contin- 
uum in AGN are well constrained by existing data, the 
soft X-ray continuum shape is often quite poorly defined 
for individual QSOs and therefore places only limited con- 
straints on BBB models. These limitations have been due 
to the low spectral resolution, signal-to- noise ratio and/or 
narrow spectral bandpass provided by previous X-ray in- 
struments, thereby making it difficult to quantify absorp- 
tion/emission features which may contaminate the X-ray 
spectrum and hence distort the derived continuum shape. 
The advent of XMM-Newton with its broad bandpass, 
high throughput and good spectral resolution now pro- 
vides the opportunity to obtain high-quality X-ray spectra 
and hence make progress in defining the intrinsic shape of 
the X-ray continuum in AGN. 

In this paper we present XMM-Newton observations 
of PKS 0558-504, a low-redshift quasar (z = 0.137; 
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Fig. 1. A Full- Window mode EPIC MOS-1 image of 
PKS 0558-504 taken during XMM-Newton orbit 42. A 
logarithmic intensity scale has been used to show the point 
spread function. No other strong X-ray sources are seen 
within the 30 arcminute field of view. 



M v « -24.5) identified during the HEAO-1 survey. It 
has quite narrow Balmer lines (FWHM = 1500 km s _1 
Remillard et al. 



1986), which has led to it being routinely 



described as a Narrow-Line Seyfert 1 (NLS1), although we 
note that the UV lines, such as Lya, are broader (FWHM 
w 3800 km s" 1 ; Kuraszkiewicz et al. pOOO) ). PKS 0558-504 
is unusual in being one of the few radio- loud NLSl-type 
galaxies (Schartel et al. 1996). We find that the X-ray 



spectrum of PKS 0558—504 is dominated by a strong soft 
X-ray excess, which may be the high-energy tail of the 
BBB. The X-ray spectrum is not contaminated by either 
cold or warm (ionised) absorption features. Soft X-ray nar- 
row emission-lines are also excluded at a high significance 
level. Thus, we claim to have detected the featureless X- 
ray continuum emission from a luminous quasar. 



2. XMM-Newton Observations 

We present observations of PKS 0558—504 taken with the 
European Photon Imaging Camera (EPIC) during orbits 
42 and 84 (1-2 March 2000 and 24-25 May 2000) of the 
calibration and performance verification phase of XMM- 
Newton. The EPIC instrument is described in Turner et 



al. Q200l| ) and Strueder et al. Q200l| ). During both orbits 
the EPIC MOS cameras were operated in several differ- 
ent observing modes. Here we include only those obser- 
vations taken with the medium filter, for which the total 
EPIC MOS exposure times for orbits 42 and 84 are 85.5 
ksecs and 53.7 ksecs respectively. Only the full- frame mode 
EPIC PN data for orbit 84 were used, for which the total 
exposure time is 10.2 ksecs. 



For operational reasons, the EPIC observations in 
each orbit were divided into several separate exposures 
of varying lengths. These exposures were reduced using 
the XMM-Newton SAS (Science Analysis Software) pack- 
age. The event lists output from the standard emchain 
and epchain scripts were further filtered for each EPIC 
exposure using the SAS xmmselect task. For the MOS 
exposures only X-ray events corresponding to patterns 0- 
12 (similar to ASCA event grades 0-4) were selected. Only 
pattern (single pixel) events were selected for the PN. 
Known hot or flickering pixels and electronic noise were 
rejected using the SAS and the low-energy cutoff was set 
to 200 eV. One of the EPIC MOS-1 images from orbit 42 
is shown in Fig. [|. 

Source spectra and light curves were extracted from 
each EPIC exposure using a circular region of diameter 1 
arcminute centred on the source position. PKS 0558—504 
is by far the brightest X-ray source in the 30 arcminute 
EPIC field-of-view. The EPIC cameras are subject to non 
X-ray background events due to particles which can vary 
significantly during an exposure. Background spectra and 
light curves for use in later data analysis were derived 
in a similar way to the source data using an offset re- 
gion close to the source. We note that the observed count- 
rate for PKS 0558—504 implies little effect due to pile-up 
and no significant spectral differences were found between 
observations taken in different MOS window modes. The 
Xspec vll.O software package was used to calibrate the 
background-subtracted EPIC spectra using the most re- 
cent camera response matrices derived from ground-based 
and in-orbit data. Before further analysis the EPIC spec- 
tra were binned to give a minimum of 20 counts per bin. 

Refle ction Grating Spectrograph (RGS; den Herder et 
al. 2001) exposures from orbit 84 were also reduced to 



search for weak/narrow soft X-ray absorption or emission- 
line features. Three exposures were analysed, two from 
RGS-1 (each of 25 ksecs) and one from RGS-2 (31 
ksecs). The RGS data were first processed using the SAS 
RGSPROC script. The merged RGS event lists were then 
fed to the SAS evselect and xmmselect tasks, and 
both source and background spectra for each exposure 
were extracted interactively. These spectra were then cali- 
brated using response matrices derived from the individual 
RGS exposures using the SAS rmfgen task. 



The Optical Monitor (OM; Mason et al. g001[ ) imaging 
data for orbits 42 and 84 were reduced using standard 
SAS tasks and calibrated with the instrument response 
information for each filter. 



3. Variability 

No substantial variability was detected within individual 
EPIC exposures on timescales less than a few minutes, 
so the data were rebinned and assembled into orbit light 
curves by combining the individual light curves for each 
exposure in order to search for longer duration events. 
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Fig. 2. The EPIC MOS-1 + MOS-2 0.2-10 keV light curve 
of PKS 0558-504 for orbit 84 in 125 second time bins. 

The PKS 0558-504 broad-band EPIC MOS light curve 
for orbit 84, binned into 125 seconds bins, is shown in 
Fig. | 

Strong broad-band X-ray flux variability is seen dur- 
ing both orbits, with similar rates of increase and decrease. 
The largest change in count rate occurs during orbit 84 - 
an increase of 35% in 3000 seconds. Based on the spectral 
fits discussed below and using H Q = 50 km s _1 Mpc~ 4 
and qo = 0.5, this corresponds to a rate of increase in 
the 0.2-10 keV luminosity of AL/At = 5 x 10 41 erg s~ 2 . 
Assuming photon diffusion through a spherical mass of 
accreting matter in which the opacity is dominated by 
Thomson scattering, the observed luminosity change im- 
plies a minimum efficiency of converting matter into en- 
ergy of rj = (AL/At)/(2 x 10 42 ) (Guilbert, Fabian & Rees 



1985 ) . In the simplest interpretation, the change in lumi- 



nosity observed during orbit 84 therefore implies rj = 0.25, 
an efficiency which exceeds the theoretical maximum for a 
non-rotating (Schwarzschild) black hole but is consistent 
with that for a maximally rotating Kerr black hole. Some 
contribution of beamed X-rays cannot be ruled out, but 
the shape of the continuum (discussed below) is not con- 
sistent with that of a typical beamed source. A tendency 
for the X-ray continuum to harden slightly when brighter 
is seen when comparing the data from orbits 42 and 84. 
The X-ray variability of PKS 0558—504 is discussed in 
more detail in a separate paper (Gliozzi et al. [2001 ). 

During a Ginga observation, the PKS 0558—504 con- 
tinuum appeared to vary by 67% in 3 minutes (implying 
■q ~ 2; Remillard et al. 1991). Given the high count-rates 



and well determined background of the XMM-Newton 
observations, we could have detected such a flare on 
timescales as short as a few tens of seconds, but none were 
observed. ROSAT HRI observa tions also failed to detect 
such rapid flaring (Gliozzi et al. 2000|) , but did show lumi- 
nosity variability rates (allowing for their use of Hq = 70 




0.5 1 2 
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Fig. 3. The EPIC MOS spectra of PKS 0558-504 for or- 
bits 42 and 84. A power-law has been fitted to the 4- 
10 keV data and extrapolated to lower energies. A broad, 
strong soft X-ray excess is clearly seen extending to ~ 3 
keV. 

km s" 1 Mpc" 1 ) which imply accretion efficiencies simi- 
lar to that derived from the XMM-Newton data. Without 
confirmation using an imaging detector it is impossible to 
confirm or deny the association of the Ginga flare with in- 
trinsic variability in PKS 0558—504, so we do not consider 
the Ginga flare any further in this paper. 

4. Spectral Analysis 

Previous analyses of ASCA observations of PKS 0558—504 
imply an X-ray spectrum which is softer at lower energies. 
Fixing the absorbing column densit y at t he Galactic value 



of 4.4 x 10 20 cm" 2 , Gliozzi et al. (g000|) find T = 2.99 ± 



0.09 over 0.1-2.4 keV from ROSAT data. ASCA da ta giv e 
T = 2.26 ± 0.03 over 0.6-10 keV (Vaughan et al. |1999[ ), 
whereas the Ginga data reveal a typical AGN power-law 
spectru m wi th T = 1.89±oq 4 over 2-18 keV (Lawson & 
Turner \L997\ ). Similarly, Comastri ( [2000D finds r = 1.97 ± 
0.04 over 2-10 keV from BeppoSAX data. 

As noted above, no large (^ 10%) changes in spectral 
shape were observed during each orbit. Thus, to improve 
the signal-to-noise ratio the MOS data were combined to 
form a mean spectrum for each orbit before detailed spec- 
tral analysis. 

The background-subtracted, summed MOS observa- 
tions provide the highest S/N data for PKS 0558-504 
to date. These spectra are very poorly fitted by a single 
power-law (x 2 = 2.75) due to the presence of a large, 
broad soft X-ray excess which dominates the X-ray emis- 
sion from PKS 0558-504 in the XMM-Newton bandpass. 
To illustrate this feature, Fig. [| shows the MOS spectra 
for each orbit with a power-law fitted to the 4-10 keV 
band, allowing for Galactic absorption (fit 1 in Table 1 - 
all the fits are summarised in Table 1). The power-law is 
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Fig. 5. The broad-band spectral energy distribution of 
Fig. 4. The multiple blackbody plus power-law fit to the PKS 0558-504 from OM and EPIC data, 
combined orbit 42 and 84 MOS data. 



fitted above 4 keV as the soft excess appears to extend to 
~ 3 keV. The PN spectrum for orbit 84 is consistent with 
the MOS allowing for the current calibration uncertain- 
ties. 

No evidence is found for a significant change in X-ray 
spectral shape between the two orbits, but there is a small 
decrease in the mean X-ray flux of w 18% from orbit 42 to 
84. The OM data show that the UV flux also decreased by 
w 10% from orbit 42 to 84. For orbit 42 the 0.2-2.0 keV 
and 2-10 keV fluxes are 2.44 x lO" 11 and 1.01 x lO" 11 
erg cm" 2 s^ 1 respectively, corresponding to luminosities 
of 3.7 x 10 45 and 0.98 x 10 45 erg s _1 . 

The soft X-ray excess in PKS 0558-504 is far too 
broad to be fitted by a single temperature blackbody. To 
parameterise the spectra with a relatively simple model, 
we fitted a combination of blackbodies and a power-law. 
At least three blackbodies are required of 80, 175 and 
475 eV, to give an acceptable fit (xl — 1.2; fit 2 in Ta- 
ble 1). Other than the relative normalisation (consistent 
with the small overall change in flux), the derived power- 
law indices and blackbody temperatures required to fit 
the mean spectra from orbits 42 and 84 are statistically 
consistent. Therefore, we combined the entire MOS data 
before proceeding. The multiple blackbody and power-law 
fit to the combined MOS data is shown in Fig. [|. The high- 
energy power-law slope in the multi-component model is 
consistent with the r = 1.9 commonly seen in radio-quiet 
AGN. 

No iron fluorescence line is detected from either neu- 
tral or ionised material. For the combined MOS data, the 
upper limits to the EW of a 6.4 keV line are < 10 eV for 
a narrow line (a = 10 eV) and < 90 eV for a broad line 
(cr = 300 eV). For a narrow ionised line at 6.7 keV, the 
upper limit is 42 eV. In addition, no evidence is found in 
the EPIC or RGS spectra for any intrinsic cold or warm 



absorption. The Ovn or Ovm edge optical depth limits 
are r < 0.1. The upper limit for an ionised iron edge is 
t < 0.2. We note that an apparent feature in Figs ||, |J 
and |^ at rs 7 keV (which is the wrong energy for an iron 
line in the QSO rest-frame) is due to noise in some of the 
MOS spectra that went to form the mean. This feature is 
not present in the PN spectra. Overall, PKS 0558—504 is 
a powerful quasar in the X-ray band, with a 0.25-10 keV 
(unabsorbed) luminosity of 4.5 x 10 45 erg s _1 mostly emit- 
ted below 3 keV. The lack of spectral features strongly im- 
plies that we are seeing the bare quasar continuum emis- 
sion. 



5. Discussion 

The multiplc-blackbody plus power-law model discussed 
above provides an acceptable fit to the EPIC spectra, but 
has no physical basis. To examine the broad-band spectral 
energy distribution (SED) of PKS 0558-504 the dered- 
dened optical and UV data from the OM (assuming Galac- 
tic E(B— V) = 0.044) together with the X-ray continuum 
derived from the combined MOS data (using the multiple 
blackbody model) are shown in Fig. |5|. 

The SED of PKS 0558-504 clearly shows an opti- 
cal/UV BBB, and the XMM-Newton data suggest this fea- 
ture extends into the X-ray band. To power an object as 
luminous as PKS 0558-504 (L bo] - 10 46 erg s^ 1 ) would 
require a black hole of mass Mbh ~ 10 8 Mq accreting at 
the Eddington rate. Based on accretion disc models (e.g. 
Czerny & Elvis 1987), the thermal spectrum of such a sys- 
tem could explain the optical/UV portion of the BBB but 
would be unable to produce the broad, soft X-ray excess 
or the high-energy power-law-like emission. 

One model suggested to explain the steep X-ray spec- 
tra of extreme soft X-ray sources, such as NLSls, is that 
these source accrete at an unusually high rate leading to a 
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Table 1. Fits to XMM-Newton data in the 0.2-10 keV band. a Powerlaw index for fit to 
4.0-10 keV. b Electron temperature of thermal source (eV). c Electron temperature of Comp- 
tonizing electrons (keV). d Optical depth. 



Rev/Camera 


Fit 


Model 


T/photon kT 


42/MOS 


1 


PL a 


2.12 + 0.15 


42/MOS 


2 


PL + 3xBB 


2.03 + 0.09 


84/MOS 


1 


PL a 


1.99 + 0.15 


84/MOS 


2 


PL + 3xBB 


1.97 + 0.12 


84/PN 


1 


PL a 


1.94 + 0.17 


84/PN 


2 


PL + 3xBB 


1.84 + 0.05 


42+84 MOS 


1 


PL a 


2.07 + 0.12 


42+84 MOS 


2 


PL + 3xBB 


1.97 + 0.10 


42+84 MOS 


3 


CompTF 


59±8 b 








59 ±% b 



electron temp (kT)/optical depth 


X 2 /dof 




241/243 


82 ± 5, 173 ± 10, 469 ± 18 


608/489 




276/258 


74 + 5, 165 + 11, 454 + 35 


622/493 




188/192 


89 ± 5, 189 ± 14, 560 ± 50 


1009/798 




319/317 


80 + 3, 175 + 6, 475 + 15 


721/556 


4.0±1j c , 2.5l?;? d 


780/559 



55t^ c , 0.7±0.4 d 



Comptonized Spectrum of PKS 0558—504 




1 2 
channel energy (keV) 
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Fig. 6. The Comptonization model fit to the combined 
orbit 42 and 84 MOS data plotted in vf„ space. 



large soft photon flux. These photons then Compton cool 
the hard photon source, generated in a hot accretion disc 
corona, resulting in a steeper X-ray spectrum (Pounds, 
Done & Osborne 1995| ). The temperature, geometry and 
covering factor of disc coronae are poorly defined from 
theory, but multiple-temperatures are likely to be required 
due to vertical temperature stratification. 

To explore Comptonization models, the results of fit- 
ting the Xspec model CompTT to the combined orbit 42 
and 84 data are given in Table 1 (fit 3) and shown in Fig. 
||. A model involving two electron populations, both scat- 
tering photons from a single temperature thermal source, 
provides a good fit to the mean X-ray spectrum. Allow- 
ing the thermal source temperature to vary independently 
does not produce a significantly better fit. The tempera- 
ture and luminosity of the soft photon source are consis- 
tent with that of a thermal accretion disc spectrum for 
a 10 8 Mq black hole accreting at close to the Edding- 
ton rate. The Comptonizing electron temperatures of 4 



and 55 keV are at the low end o f the expected temper- 
ature range (Haardt & Maraschi 1993 ), possibly due to 
the large soft photon flux Compton cooling the corona. 
We note that the exact values of the parameters for the 
two electron populations in the comptt model are quite 
poorly constrained. This is because increasing the optical 
depth slightly can harden the spectrum in a similar way 
to adding a higher temperature electron population, thus 
giving a variety of two-component models that adequately 
fit the data. 

We note that, as suggested by the power-law and mul- 
tiple blackbody fits, the high-energy tail of the Comp- 
tonized spectrum has a shape consistent with that of 'nor- 
mal' broad-line AGN power-laws (r ~ 1.9), and is signif- 
icantly harder than that derived from the ASCA data. 
These differences are presumably due to the inability of 
the ASCA data to delineate the tail of the soft X-ray ex- 
cess which extends to around 3 keV in this source. The ob- 
served mean slope of the hard power-law in AGN (r = 1.9) 
is predicted from models in which feedback occurs between 
soft photons produced by hard X-ray heating Compton 
cool a moderate optical-depth corona (Haardt & Maraschi 
1993 ). Such a model cannot explain the bulk of the soft 
X-ray luminosity in PKS 0558—504 as it dominates the 
luminosity, but may produce the residual hard power-law 
emission. 

The reduced chi-squared value for the Comptonized 
model (xi = 1-3) shows that the Comptonization model 
fit leaves some spectral curvature, mainly above a few 
keV. Some of this imperfection is due to remaining calibra- 
tion uncertainties around features such as detector Gold 
edge, but most of it is due to the Comptonization model 
simply being too simplistic to adequately fit such high- 
quality data. Some of the remaining spectral curvature 
can be reduced by including an additional Comptoniz- 
ing component (i.e. hotter electrons) but this component 
does not formally improve the fit in terms of the reduced 
chi-squared presumably because the XMM-Newton high- 
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energy cutoff is too low to adequately constrain a very-hot 
electron population. 

An alternative explanation for a high-energy 'power- 
law-like' component is a small amount of reflection from an 
ionised disc. The disc is likely to be ionised given the lumi- 
nosity of the source. Indeed, the absence of any strong iron 
Ka emission line or iron edge at high energies and the lack 
of absorption features due to moderately ionised material 
at low energies suggests that any X-ray heated disc will 
have a substantial 'skin' which is fully ionised and effec- 
tively acting as a 'mirror'. This has been suggested for high 
accretion rate sources in the context of a magnetic-flare 
model for X-ray variability (Nayakshin 200C| ). The small 
iron Ka emission line and iron edge strengths predicted 
by the magnetic-flare model, assuming r ~ 2, are consis- 
tent with the upper limits from the XMM-Newton spectra 
(Nayakshin, private communication) . What is less clear is 
if a fully ionised corona can be maintained in a source like 
PKS 0558—504 with a large soft photon flux if those pho- 
tons are generated close to the magnetic-flares. Clearly the 
XMM-Newton data for PKS 0558-504 demand compari- 
son with more detailed accretion disc models. 



6. Conclusions 

The XMM-Newton X-ray spectrum of PKS 0558-504 is 
dominated by a large, variable soft excess. The fastest ob- 
served variability implies a minimum efficiency of convert- 
ing matter to energy of 77 = 0.25, suggesting accretion 
onto a Kerr black hole. Combining the X-ray data with 
simultaneous OM optical and UV data reveal a large big 
blue bump; too broad to be due purely to a thermal ac- 
cretion disc spectrum. The X-ray spectrum can be fitted 
by a Comptonization model in which thermal photons, 
presumably from the disc, are scattered by a hot corona. 
At least two electron temperatures are required. Overall, 
PKS 0558—504 can be explained by an accretion disc run- 
ning at a high accretion rate producing copious amounts 
of soft EUV photons which are then Compton cooling a 
hot disc corona. 
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